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Abstract 



> 
in 

O ■ Realistic predictions are made for the rates of lepton flavour violating Higgs boson 

2 ' decays, r /i7, ^ 67, yU/i, 5° r/i and r via a top-down 



analysis of the Minimal Supersymmetric Standard Model(MSSM) constrained by 
O I SU (5) unification with right-handed Neutrinos and large tan /5. The third family 

'-^ ! neutrino Yukawa coupling is chosen to be of order 1, in this way our model bares 

^ . a significant resemblance to supersymmetric 50(10). In this framework the large 

^ ■ PMNS mixings result in potentially large lepton flavour violation. Our analysis 

predicts r ^ fi'y and fi ^ c-f rates in the region (10~^ — 10~^) and (10~^^ — 
10~^^) respectively. We also show that the rates for lepton flavour violating Higgs 



^ I decays can be as large as 10"''. The non-decoupling nature of — > r/i is observed 

which leads to its decay rate becoming comparable to that for r for large 

values of mo and M1/2. We also find that the present bound on Bs fifi is an 
important constraint on the rate of lepton flavour violating Higgs decays. The 
recently measured Bg-Bg mixing parameter AM^ is also investigated. 
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1 Introduction 



Over the past few years the phenomenon of neutrino oscillations and neutrino masses have 
been confirmed by observations made by the Super-Kamiokande[l] , SN0[2], K2K[3] and 
KamLAND[4] experiments. At the present time the experimental determination of oscillation 
parameters has become increasingly accurate. The resulting measurements have indicated 
large, almost bimaximal mixing in the lepton sector, hence Neutrino mixing appears to be 
markedly different from that of the quark sector of the Standard Model(SM). A combined 
Global fit to the three-neutrino framework leads to the following 2a oscillation parameter 
fit[5], 

5m?2 = 7.92(1 ± 0.09) X 10"^ eV^ at ± 2(7 (1) 
sin^^is = 0.314(11°;}^) at ±2a (2) 

Am^3 = 2.4(ll°:2^) X 10-^eV^ at ± 2(7 (3) 
sin^e23 = 0.44(ll|ji^) at ±2(7 (4) 

sin^ ^13 = 0.9lg| X 10-^ at ± 2a. (5) 

Neutrino oscillations are the first example of physics beyond the Standard Model. In 
the Standard Model of particle physics the lepton number of each generation is a conserved 
quantity. The observation of neutrino masses also marks the first evidence for fiavour viola- 
tion in the lepton sector and gives the possibihty of lepton flavour violation (LFV) among the 
charged leptons. In the supersymmetric seesaw model, flavour violation at the high energy 
scale may feed into the charged lepton sector via radiativcly induced off-diagonal elements of 
the left-slepton scalar mass, m|. Hence the large mixings observed in neutrino oscillations 
could result in large lepton flavour violating rates for such decays as r — > — > 67 and 

Supersymmetric theories contain a number of possible sources of lepton flavour violation. 
As a result, the bounds on rates of LFV are particularly restricting upon the SUSY parameter 
space. Even in the case of minimal supergravity(mSUGRA), where soft SUSY breaking 
masses and trilinear couplings are flavour diagonal at the GUT scale, the presence of right- 
handed neutrinos and the PMNS mixings are enough to radiatively induce large LFV rates. 
It is these LFV rates, induced by renormalisation group(RG) running that we wish to study 
in the present work. 

The processes — > /j7 represent the most stringent bounds on lepton flavour violating 
interactions. The present experimental bounds are as follows, 

(6) 
(7) 
(8) 

In the near future these decays will be studied at the PSI-MEG experiment and at the 
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B-factories, which will set bounds in the region of, 
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10-« [7] (10) 
10-« [7]. (11) 

When studying lepton flavour violating decay rates it is important to ensure that the rates 
for these rare lepton decays don't exceed the experimental bounds listed above. In order to 
make realistic predictions for LFV rates within a complete SUSY model, it is useful to require 
a good fit to electroweak data including, fermion masses and mixings, correct electroweak 
symmetry breaking and the correct rate for b 57. 

Recently a number of Higgs mediated lepton flavour violating processes have also be- 
come interesting. These include, -B^^ TAi[9], r — >• /x/i/i[10], r and lepton flavour 
violating decays of neutral Higgs bosons[12]. In the region of parameter space for which 
the neutral Higgs boson masses are hght, these Higgs mediated decays can possibly become 
very interesting. Enhancement by large tan adds to the possibilities even further. When 
discussing such rare decays it is very important to take into account the experimental bounds 
from li — > Ij 7, only then can meaningful predictions for such decay rates be made. Theoreti- 
cal studies of such decays with supersymmetric models have reported interesting predictions 
with rates as high as, 

BR(5° ^ r/x) ~ 4 X 10"^ [9] (12) 
BR(r iififi) ~ 10"^ [10] (13) 
BR(r ^ firi) ~ 10 x BR(r ^ /x/x/x) [11] (14) 
BR(/ ^ r/x) ~ 10"^ [12], (15) 

where 0° represents the neutral MSSM Higgs bosons, (f)^ — A^, H^, hP. 

The most exciting example of Higgs mediated decays is that of — > ////[13] in the 
MSSM. This decay is proportional to tan^ /3 and so is greatly enhanced by large values of 
tan j3. If the Higgs boson is very light then this rate could become very exciting. The present 
experimental bounds for these rare decays are the following, 

BR{Bs nii) < 0.8 X 10-^ at 90% C.L. [14] (16) 
BR{Ba nn) < 2.3 x 10"^ at 90% C.L. [14]. (17) 

Theoretically the ratio Vts/Vtd ~ 10 means that the Bg decay is more restrictive than that of 
Bd- These decays are of particular interest as the Standard model contributions enter at the 
1-loop level which makes this a good place to search for physics beyond the standard model. 

In supersymmetric models it is expected that the dominant contribution to Bg^d — ^ l^l^ could 
come from a penguin diagram mediated by Higgs bosons. In such a scenario the present 
bound on Bg —>■ will constrain the Higgs masses and in turn constrain the decay rate for 

4P ^ Til. 

We propose to study the rates of lepton flavour violating decays of r, B^^^ and MSSM 
Higgs bosons in a supersymmetric model based on SU (5) unification with seesaw neutrinos 
and large tan /3. The quark and charged lepton sectors of this theory are well constrained 
with 9 parameters determining 12 low-energy observables. On the other hand the neutrino 
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sector has more parameters than experimental observables. Therefore we must attempt to use 
bounds on lepton flavour violating decays to further constrain the parameters of the neutrino 
sector. Prom this theory we hope to make general predictions for some of the interesting 
lepton flavour violating decays mentioned above and also to study their correlations. The 
recently measured Bg-Bg mixing parameter AMg is also investigated. 

The present work differs from previous studies of lepton flavour violating Higgs decays 
in that our work uses a top-down analysis of a Grand Unified SUSY model. Also we are 
studying lepton flavour violating Higgs decays in conjunction with Higgs mediated decays, 
such as Bg —>■ for the first time. It is shown that the experimental bound on Bg — > //// 
provides important information regarding the allowed rates for lepton fiavour violating Higgs 
decays. It is also shown that the Higgs contribution to AM^ via the often neglected operator 
{bRSL){pRSL) may be comparable to the contribution from the operator {hRSL){hLSR}. 

The rest of the paper is arranged as follows. A basic introduction to the supersymmetric 
SU{5) model studied in the present work is given in section 2. A discussion of the relevant 
fiavour violating phenomenology analysed in this paper is given in section 3 with an outline 
of our numerical procedure in section 4. In section 5 we present our results and discuss their 
implications. Finally, section 6 concludes the present work. 

2 Minimal Supersymmetric SU{b) with right-handed 
neutrinos 

In the present work we shall study the MSSM+i/ij constrained at the GUT scale by SU (5) 
unification. Let us first outhne the general features of this model and then proceed to give 
details of the simplifying assumptions made during this analysis. In Minimal SU{5) the 
matter superfields of the MSSM are contained within the representations; 10 = {Q, If^, E'^), 
and 5 = (L, D'^). The inclusion of right-handed neutrinos into the theory requires a further 
singlet superfield, 1 = {N'^)- In addition there are 5 and 5 Higgs representations, H and H. 
These Higgs representations contain the MSSM Higgs doublets H^. 

Using these matter superfields we can construct the Yukawa section of the superpotential 
as follows. 

The first term of eq. (18) gives rise to a symmetric up-quark Yukawa coupling Yu and the 
second term provides both down-quark and charged lepton Yukawa couplings. As a result 
we have the GUT scale relation, — , hence there exists GUT scale h — t unification^. 

The final two terms in the superpotential of eq. (18) are responsible for the neutrino 
Yukawa couphng, Yi, and the right-handed neutrino Majorana mass, Mr. These final two 
terms combine to produce the light neutrino mass matrix via the type I seesaw mechanism, 

niLL= -vIY^Mr-'yJ. (19) 

^Clearly this relation cannot hold for the first two generations of down quarks and charged leptons. In such 
a case it is possible to invoke corrections due to non-renormalisable operators or suppressed renormalisable 
operators[15]. 
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Our next step is to rotate away the non-physical degrees of freedom. We begin by rotating 
the 10 and 5 representations such that is diagonal. As a result the coupling F„ is solely 
responsible for the CKM mixings. Therefore, from eq. (18) we can see that Yy, is diagonahsed 
by rotating the 10-dimcnsional representation by the CKM matrix. The Icpton sector is more 
complicated due to the structure of the see-saw mass matrix in eq. 19. Rotating the singlet 
neutrino such that the combination, Yj}Y^, is diagonal we can write. 



Yu = V^^^Y^'^^VcKM 
— 

Y^ = UYf''^ (20) 
Mr = WM^^^W^. 

Here U and W both represent unitary rotations. In the case were, W — 1, the mixing matrix 
U is identified as the PMNS matrix^. For simphcity we assume that W — 1 and that all 
parameters are real. ^ 

In the neutrino sector we have 9 parameters at the GUT scale, but only 4 low-energy 
observables. The bound on the undetermined neutrino mixing angle, ^13, also constrains the 
GUT scale mixing to be small. As a result and for simplicity we shall assume ^i3(Mgut) = 
in our present study. Throughout the analysis we have made the assumption that both 
the neutrino Yukawa and Majorana matrices have a hierarchical form with Y^^^^^ ~ 1 and 
Mr^ ~ 10^^ GeV, as shown in eq. (21) and (22). We choose such a hierarchy in the neutrino 
Yukawa sector with one eye on b — t — t Yukawa unification and SO{10) with tan /3 ~ 50. 

ydiag ^ ydiag ^ ydiag ^21) 

M^fs < M^;"s < Mjf ~ 1014 GeV . (22) 

This type of hierarchy naturally leads a light neutrino scale ~ 0.1 eV. It seems that we are 
still left with a great deal of freedom in the neutrino sector with four parameters at the GUT 
scale determining only two mass squared differences at the weak scale. It is worth noting that 
if we had studied SO {10) unification then the additional symmetry would provide a relation 
between each of the fermionic Yukawa couplings and so the neutrino sector would then be 
more predictive. Here we study SU{5) as it is the prototypical unified group and in many 
ways it is the minimal unified theory. In SU{5) the freedom in the neutrino sector must 
be constrained by bounds on the rates of lepton flavour violating decays such as — > ej, 
r — > //7 and r — > 67. Prom eq. (20) the neutrino Yukawa couphng at the GUT scale, 
assuming bimaximal neutrino mixing and ^13 = 0, takes the form. 



Y 



/ _i_ "V^diag _i_ -y^diag n 

_ 1 -i^diag 1 Y'diag J_ -y^diag I /oon 

1 y^diag -1 "y^diag _l_ y diag 

\ 2 1^1 2 ^2 a/2 ^^3 



Eq. (23) shows that the size of the off-diagonal terms in the neutrino Yukawa matrix 
are determined by the diagonal elements, Y^^^^, Yv2^^ Y^z^^- Such terms will induce 

^Such a case arises naturally in U{1)f flavour models where the U{1)f charges of the right-handed 

neutrino is responsible for the right mixing in Yi, as well as the mixing in Adn. 

•^A similar SU{5) framework has previously been used to study fi ^ and hadronic EDMs etc. [16]. 
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off-diagonal elements in the slepton mass squared matrix with the approximate form, 



(A,4), ~ _ h(iVWM«) ^^^^^g ^ 2AI){KYJ),, r + j. (24) 



Here mo and are the usual universal soft scalar mass and trilinear couplings commonly 
present in minimal supergravity models. Finally this leads to radiatively induced lepton 
flavour violation through Feynman diagrams with sleptons inside the loop [17]. Clearly 
eq. (24) tells us that it is the size of Y^l^'^ that will affect the rates for — > l^^. For example 
with, y,^3 > y,^2, it will be Y^^ which contributes strongly to the rate of r — > //7. Hence the 
bounds on ^ 67, r ^ ii'^ and r 67 will act as a constraint on the size of Y^^^^. During 
our study we have chosen a neutrino Yukawa hierarchy of the approximate form, 

y diag . ydiag . y diag _ : : 1 (25) 

where A = 0.22. This neutrino Yukawa hierarchy is similar to the hierarchy in the up-quark 
Yukawa coupling, which naturally arises in simple 5*0(10) models [18]. With, Y^^ ~ 1, it is 
possible that /x-r lepton flavour violating rates will be near their experimental limits. This 
allows us to analyse predictions for other LFV decays while the rates for r — > /X7 are at the 
maximally allowed levels from experiment. In this way our predictions will be optimistic 
while remaining realistic. 

Let us now examine the soft SUSY breaking section of the Lagrangian. For minimal 
SU (5) the relevant section of the soft supersymmetry breaking Lagrangian takes the following 
form, 

-Cson = ^10i(mio)iiloJ- + §- (m5)ij§j + ii(m^)ij it + m|h^h + m| 

+ ^ fOi {Au)ij 10,- h + {Aa)ij 10,- h + (^jv)i,- h h + ^^24 2424. (26) 
o z 

Here 10, 5, 1, h and h are the scalar components of the fields, 10, 5, 1, H and H. The final 
term of eq.(26) is the soft mass for the SU{5) gaugino 24, which is the scalar component of 
the gauge boson adjoint representation, 24. Below the GUT scale the soft SUSY breaking 
terms of the effective theory are given by. 



soft 



+ QMd)ijdjhd + Li{A^)ijejhd + Li{Ar,)ijnjh^ (27) 

where, Q, L, u, d, e, n, hu and hd are scalar components of the superfields Q, L, W^, D'^, E^, 
N"^, Hu and Hd respectively. B, W and g are the gauginos of the MSSM contained within 
the adjoint 24. Here the gaugino soft masses are related by the SU{5) symmetry such that. 
Ml = M2 = M3 = M24. From eq. (27) we see that tjiIq provides the soft scalar masses, uiq, 
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ml, ml, and mj gives the masses ml, and ml. Hence we have the following conditions, 

TTT'Q = ^lo^, = ^?0' = 

(28) 

ml = m|, m|; = m^Q. 

In the present work we are interested in the renormalisation group effects of off-diagonal 
elements of the Yukawa couplings. Due to the large mixings in the PMNS matrix, these off- 
diagonal elements can be large. We shall also assume the GUT scale flavour blind boundary 
conditions, 

{mlo)ij = {ml)^J = (rr4),j- = 

= AoY:^, = A^^, A'J = AoV:^, Al= = A^Y^^, (29) 
Ml = M2 = M3 = Ml/2, 

where, mg, Aq and M1/2 are the universal soft scalar mass, trilinear coupling and soft gaugino 
mass respectively. Throughout our analysis we shall also assume = at the GUT scale. 
Under these assumptions the soft Lagrangian becomes flavour blind and hence contains no 
new lepton flavour violating sources. Therefore, in the soft SUSY breaking Lagrangian we 
have just five parameters, mo, Afi/2, ^0, nT'l and m'g. The assumption of universality ensures 
that flavour violation induced by off-diagonal scalar masses does not become too large. Also it 
will allow us to assess the renormalisation group (RG) running effects of off-diagonal Yukawa 
couplings. 

The SU (5) superpotential, as well as providing the Yukawa couplings of the MSSM+z//^ of 
eq. (18), also includes quark and lepton couplings to the coloured Higgs. Such couphngs may 
also induce flavour violating interactions in the squark and slepton mass squared matrices 
[16]. Here we shall ignore these couplings as they may also induce Proton decay at an 
unacceptable rate. Therefore, it seems to be an acceptable assumption for such couplings to 
be suppressed. 

The next section is devoted to a discussion of the flavour changing neutral current and 
lepton flavour violating phenomenology studied in our analysis. 

3 LFV and FCNC Phenomenology 

Observations of neutrino oscillations at SuperK [1], SNO [2], K2K[3] and KamLAND [4] 
imply the existence of massive neutrinos with large solar and atmospheric mixing angles. The 
small neutrino masses are most naturally explained via the seesaw mechanism with heavy 
singlet neutrinos. Even in a basis where both and Mr are diagonal in flavour space Yi, is 
always left as a possible source of flavour violation. In SUSY models this flavour violation can 
be communicated to the slepton sector through renormalisation group running. The initial 
communication is from running between the GUT scale and the scale of M^. Although 
the scale Mji is far above the electro-weak scale its effects leave a lasting impression on the 
mass squared matrices of the sleptons. Subsequently flavour violation can enter into the 
charged lepton sector through loop diagrams involving the sleptons and indeed such effects 
have been used to predict large branching ratios for r — > //7 and // — > 67 within the MSSM 
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[19, 20, 21, 22]. 

In the Standard Model Flavour changing neutral currents are absent at tree-level and 
only enter at 1-loop order. As a result FCNCs are heavily suppressed. In extensions of the 
SM, the MSSM for example, there also exist additional sources of FCNC. A clear example 
comes from the mixings present in the squark sector of the MSSM. These mixings will also 
contribute to FCNCs at the 1-loop level and could even be larger than their SM counterparts. 
An example that we shall study in this work are the flavour changing couplings of neutral 
Higgs bosons and the neutral Higgs penguin contribution to such decays as 5° ^fi. When 
examining such flavour changing in the quark sector it is obviously important to consider 
the measurement of the decay rate for 6 — 57. 



3.1 Higgs mediated Flavour Changing Neutral Currents 

It has recently been pointed out that Higgs mediated Flavour Changing Neutral Cur- 
rent(FCNC) processes could well be among the first signals of supersymmetry[13]. In the 
MSSM radiatively induced couplings between the up Higgs, Hu, and down-type quarks may 
result in flavour changing Higgs couplings. In turn this will lead to large FCNC decay rates 
for such decays as Bg fi^fi~. In the Standard Model the predicted branching ratio for 
Bg is of the order of 10~^, but in the MSSM such a decay is enhanced by large tan /? 

and may reach far greater rates. It is also possible to extend this picture to the charged 
lepton sector where similar lepton flavour violating Higgs couplings are radiatively induced. 
These LFV couplings can then result in such Higgs mediated decays as r — /i/i/i, the neu- 
tral Higgs decays (p^ Tfi and even combined with the afore mentioned FCNC couplings 
resulting in the decay Bg — » fir. 
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Figure 1: Examples of Feynman diagrams contributing to the effective h'^sH^ coupling. 
Black dots indicate flavour-changing vertices while crosses stand for mass insertions for 
interaction eigenstates. 



In this paper we shall be studying these decays in the framework of supersymmetric 
SU{5) where all the tree-level MSSM couplings are determined at the high energy scale. 
In order to make realistic predictions for 23 flavour transitions it is important to analyse 
a complete theory such as this one. In this way we are able to constrain the 23 sectors of 
the theory by demanding the theory reproduces correct quark and lepton masses, CKM and 
PMNS mixings, as well as rare decay rates such as 6 — 57 and r fi'j. Only through 
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these strict constraints can realistic predictions be made for the potentially interesting decay 
channels discussed earlier. Let us now review the sources of these flavour changing Higgs 
couplings. 

As we mentioned above, loop diagrams such as those shown in fig. 1 induce flavour 
changing couplings of the kind, h'^sH'^* . Wc can write such effective couplings as functions 

fd and Qd in the Lagrangian of eq. (30). Here, fd and Qd are matrices in flavour space. The 
function Qd represents the effective coupling of the up Higgs resulting from diagrams such as 
those in fig. 1. The other coupling, /d, is the corresponding Hd effective vertex induced by 
similar diagrams with Hd as one of the external legs. 



-1 



:(0) 

R 



(30) 



In eq. (30) c?^''^ and Y^'''''^"'^ are the tree-level mass eigenstates and Yukawa coupling of the 
down-type quarks. Going beyond tree-level introduces the effective vertices, fd and Qd, which 
therefore also provide additional mass term contributions as follows. 



= -d 



(0) 
R 



m 



(O)diagt 



+ fd'"d + gW 



d 



(0) 



(31) 



If Vu ^ Vd this will lead to sizeable corrections to the mass eigenvalues [23] and mixing 
matrices [24]. Furthermore the 3-point functions in eq. (30) and mass matrix in eq. (31) 
can no-longer be simultaneously diagonalised [25]. Hence, beyond tree- level we shall have 
non-diagonal Higgs couplings in the mass eigenstate basis. 

We can write eq. (30) as. 



-dR 



^(O)diagt , ft , „t!^ 
'^Jd^yd^, 

Vd 



Vd 



9'd K 



d 



(0) 



(32) 



Notice that the first bracket of eq. (32) is in a form similar to that of the mass matrix 
and therefore is diagonal when d^l^\ arc rotated into the 1-loop mass eigenstates, d^^"^^ = 



L,R(1) ,(0) 



dj^ j^. The second bracket on the other hand is not in a similar form and so will be 



non-diagonal in the 1-loop mass eigenstate basis. This second bracket therefore becomes a 
flavour changing Higgs coupling. 



■'FCNC — 



-d 



(1) 

Ri 



Vd^'''9\[HT-':^H^d 



dl] + h.c. 



(33) 



It is clear that the origin of this flavour violating source is the mismatch between the 1- 
loop mass matrix and the couplings to the Higgs and Hd- From eq. (33) we see that 
this mismatch is caused by the interaction dp^H^*d^i^^ not present at tree level. Notice that 
the Cfcnc vanishes if Qd = 0. Eq. (33) also shows that these flavour changing couplings 
are enhanced by an explicit factor of tan/3 on top of any tan/3 scaling present in Qd- In 
the leading order in tan/3 the Qd matrix can in fact be related in a simple way to the 
finite non-logarithmic mass matrix corrections, {gd)ij = {^tTid^"'^^'^)ij/vu, computed for the 
first time in [24]. Decomposing the Higgs fields in terms of the physical Higgs boson states, 
i/o = Vu + {H^Sa + h^Ca+iA^cp+iG%(3)/^/2 and = va + iH^Ca-h^Sa+iA^sp-tG^cp)/ V2 
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we can write, 

- -H', = ^ - [H's^., + h'c^_, - ^ , (34) 

where Sa = sin a, Ca = cos a, etc. We can thus identify effective vertices IirSlH^, i>RSLh^ 
and IirSlA^ involving 6 to s transitions mediated by neutral physical Higgs states. We note 
that with large tan (3 the coupling to the pseudoscalar is always large while the CP-even 
states, and have couplings which depend on the CP-even Higgs mixing angle a. The 
Goldstone mode is cancelled in the equation above and thus the effective vertex with the Z 
boson is absent at this level. 




Figure 2: Higgs penguin contribution to the flavour changing neutral current process, 
Bs — >■ fi~^fi^. The coupling gi,s is an effective vertex generated from loops where the heavy 
SUSY partners have been integrated out. The mediating stands for neutral Higgs mass 
eigenstates, A^. 



In the MSSM with large tan f3 it is possible for the dominant contribution to Bg £^i~ 
to come from the penguin diagram where the dilepton pair is produced from a virtual Higgs 
state [13], as shown in fig. 2. After the SUSY partners are integrated out we are left with 
the effective vertices determined above. Thus in combination with the standard tree-level 
term Cuh = —ydR^^bH^ + h.c. the dominant tan/3 enhanced contribution to the branching 
ratio turns out to be 
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where the matrix 5m|^ is in the {d^l^^^} basis, and is defined by 



M- 



170 



tan/? 
50 



(35) 



(36) 



nib is the b quark mass at scale Mz in the effective SU{3)c x U{l)em theory, the constants 
are Mioo = 100 GeV and Muo = 170 GeV. The numerical factor in eq. (35) arises from: 



1.75 X 10" 



tbPbMI 0.042 0.0311^50^ 



1287r 



(37) 



In the large tan/? limit the Higgs Double Penguin(DP) contribution to B^-B^ mixing is 



10 



dominant [26]. For convenience let us write the FCNC coupling of eq. (33) as, 



Cfcnc 
Here we have written, 



X 



5° 
RL 



X 



5° 
LR 



dfij S . 



X 



RL 




(38) 



(39) 



where, A 



so 



Sa-(3, Ca-f}, -«), for = (if°, A°) . It is clear that the FCNC couplings 



are related as, [Xrl] = [Xlr] ■ In general we should also notice that, [Xrl] 



Hence, in the case of, (z, j) = (6, s), we have [Xrl] 



bs 



J_ 



RL 



LRlbs- 



Sr{sl) 



^LR{RL)ll>a 




■^LR{RL)ll>a 



bdbR) 



(a) 




Figure 3: Higgs Penguin contributions to Bg-Bs mixing where, (a) contributes to the operator 
Qi^^iQi^^) and (6) contributes to the operator Q^^. 



Following this notation we can write the neutral Higgs contribution, fig. 3, to the AB = 2 
effective Hamiltonian as. 



njAB=2 
^eff 



1 [X'RLlbsjX'^ljbs ^SLL I 1 \XLR\hs [^fj?,]fes ^SRR 

2^ -Ml ^1 +2^ -M| 

_,Sr^ [XRi\bs[Xi^bs ^LR (An.\ 



where we have defined the operators. 



Qf^ = (bPLs){bPLs) (41) 
Qf«« = {bPns)(bPRs) (42) 
Q^^ = ibPLs)ibPRs) (43) 



It is common at this point to notice that the contribution to Q2 is dominant over 

due to a suppression from the sum over Higgs fields, T~ = i -j^ + — -gr j . The 

(^2 ^2 \ 

-jir- + -j^ + jfTj ■ It turns out that this 

typically results in a suppression factor, ^ < JF^/JF+ < Recalling that [X^^ijjbs ~ 
■^[XRijbs, it may be possible for the Qf^^ contribution to give a significant effect. On the 
other hand, the contribution to Qf^^ is heavily suppressed. 
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Following the above conventions we can write the double penguin contribution as, 

= -^MbJbA 1^ (44) 

s 

Si \vS 



■^MsJbA 1^ Ml 



5 

Here P^^^ — —1.06 and = 2.56, include NLO QCD renormalisation group factors [26]. 
After taking into account the relative values of J^^, the two P's and the factor of 2 in eq. (44), 
we can see that there is a relative suppression, ^ < AMf^/AMf-^ < |. As a result the 
contribution M^^ from the operator Qf^^ may not be negligible. This relative suppression 
shall be discussed further in the following section. 

Unfortunately there is presently a large non-pcrturbative uncertainty in the determina- 
tion of Jb^- Two recent lattice determinations provide [27, 28], 

^ 230±30MeV (45) 
/B°f = 259±32MeV, (46) 

which in turn give different direct Standard Model predictions for AMf^, 

^^SM'o4 ^ 17.8 ±8ps-^ (47) 
^^SM'o5 ^ I9.8±5.5ps-^ (48) 

The recent precise Tevatron measurement of AMg is consistent with these direct SM predic- 
tions but with a lower central value [29] , 

AMf = n.Sltoil ± 0.07ps-^ (49) 



3.2 Higgs mediated lepton flavour violation 

We have seen that flavour changing can appear in the couplings of the neutral Higgs bosons 
and is enhanced by large tan/3. In the quark sector interactions of the form, dudLH^* , are 
generated at one-loop [23, 24] and at large tan /3 can become comparable to the tree-level 
interaction, dRdiH^- These two contributions cannot be simultaneously diagonalised and 
lead to potentially large Higgs-mediated flavour changing processes such as Bs — > [13]. 
Similar Higgs-mediated flavour violation can also occur in the lepton sector of SUSY seesaw 
models through interactions of the form crClH^*- This leads to the possibility of large 
branching ratios for Higgs-mediated LFV processes such as — > r//, r — > 3// and lepton 
flavour violating Higgs decays [9, 10, 12]. 

In the SUSY seesaw model the neutrino Yukawa coupling induces flavour violation in the 
slepton sector. This in turn results in flavour mixing among the charged leptons via loops 
including sleptons. As we did with down-quarks, we can write effective couplings fe and Qe 
after heavy sparticles are integrated out of the Lagrangian so that. 
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Here, e^''^, represent the tree-level mass eigenstates. There are also mass term contributions 
of the same form but with the Higgs fields replaced by their VEVs. 



-X- 



/ ~N - / \ 

/ \ / \ 

' X ^ — ' , X _ — X — ' — 

Hi' 

/ \ / \ 

/ \ / \ 

Tr < ^i-L Tr ' ^ ti-L 

' X — ^ — X — ' — ' — X — ^ — X — ' — 

Hi Hu\ W Hd H^\ B 

Figure 4: Diagrams that contribute to the coupling tr^lH*. Black dots indicate flavour- 
changing vertices while crosses stand for mass insertions for interaction eigenstates. 



Once again we see that at tree-level fe = 9e = ^ and the Yukawa couplings and mass 
matrix for the charged leptons can be simultaneously diagonalised. At one-loop level /e and 
Qe are to be computed and it follows that the one-loop 3-pt couplings and mass matrices are 
no longer simultaneously diagonalisable. The cause of this is the term, e^^glHl*ef\ Fig. 4 
shows the dominant contributions to the effective vertex g^.. 



The flavour changing part of the Lagrangian therefore takes the form. 



with. 



LFV 



0* 



Vd 



V, 



e(l)t 



(51) 



Here the matrices Vf'^^ rotate the fields from tree-level mass eigenstates e^l''n to the one-loop 



Vd 



1 1 

^2 cp 



(52) 



(0) 



L,R 



L,R 



mass eigenstates, e^^^^^ = Vljl ^l^r- As we did for the down quarks we can relate the matrix 
Qe, at leading order in tan/5, to the finite non-logarithmic corrections to the charged lepton 
mass matrix, {ge)ij = {6mf''^'')ij/vu. 

Possibly the most interesting application of such lepton flavour violating couplings is in 
the decays of MSSM Higgs bosons. The measurement of these decays at futures colliders 
could allow the direct measurement of the Higgs LFV coupling. The lepton flavour violating 
Higgs couplings discussed above facilitate such flavour violating Higgs decays. Using the 
notation used earlier we can write the partial widths of the lepton flavour violating Higgs 
boson decays within the MSSM as. 



r 



1 {Sm,)l 2 



1677 



Ms (l-Xi- Xj) ./(I - (xi + Xj)^)il - (xi - x,)2) (53) 



[sm{a- /3), cos{a- (3), -i] for ^ = [ A° ] . 
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Here, = h^, A'^ represent the three physical Higgs states, Ms = Mh, Mh, Ma are 
their masses, U = t, fi,e are the three charged lepton states and Xi = {mijMsY- In order to 
compute the branching ratio for these decays we must also compute the full width of each 
of the MSSM Higgs bosons. Throughout our analysis we shall be using the HDECAY[30] 
numerical code for the calculation of the full widths. A few authors [12] have already made 
predictions for lepton flavour violating decays within the MSSM-seesaw model. Our work 
differs from those as we are using a top-down approach including a full analysis of electro- 
weak data. We are also studying 0° r/x and Bg — > /x/x together for the first time. 



Figure 5: Lepton flavour violating decay of MSSM Higgs bosons via the effective coupling 



We can also make use of the LFV Lagrangian term of eq. (51) to study the Higgs mediated 
contributions to the process r — >■ /x/x/i for example. The dominant Higgs contributions will 
come from the penguin diagram shown in fig. 6a. There is of course a contribution from the 
photon penguin which can be related to the branching ratio of r ^ fi'-f as [17], 



Br(r ^ i^-f) 



0.003. 



(54) 



This relation is model independent and so it is possible for us to apply a rough bound 
of Br(r 3/i)^ < 9.3 x 10~^° on the photon penguin contribution from the present bound 
of Br(r fi'-f) < 3.1 X 10~^ [7]. The present experimental bound of Br(r —>■ 3fi) < 
1.9 X 10^'' [31] is at present a few orders of magnitude from this level, but in the future any 
measurements that significantly deviate from this relation, eq. (54), would be a clear signal 
that additional contributions, such as Higgs mediation, are present. Hence predictions for 
the Higgs mediated contribution Br(r —>■ ^^n) > 10"^ would indicate such a possibility. 

In the MSSM with large tan (3 the dominant contribution to the branching ratio of r ^ 
/i/i/i turns out to be. 



Br(r ^ 3/i) 



40967r3 



(^me)23 



Ml, 



■ (55) 



Here is the lifetime of the tau lepton and is the Yukawa coupling of the muon. The 
present experimental bound for this decay is as follows. 



BR(r 3/i) < 1.9 X 10"^ at 90% C.L. [31] 



(56) 
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It has also been noted [11] that the related process r — > /ir^ would in fact yield a larger 
branching ratio. The enhancement of this process comes from a factor of 3 for colour and a 
factor of {nis/Tn^Y the Yukawa coupling. In addition to this the crossed diagram in the 
muon case lowers the rate by a factor 3/2, hence the overall enhancement is by |^ ~ 10. 
The present experimental bound is given as, 

BR(r iiT]) < 3.4 X 10"^ at 90% C.L. [32]. (57) 

and so it is clearly more constraining than t ^ 3fi. This bound implies a related constraint 
on the T —>■ 3fi decay of the order, BR(r — »• 3/x) < 3.4 x 10~®, which is now only about an 
order of magnitude above the estimated photon penguin limit. It seems that there may not 
be much room for any new physics to show itself in this rare decay. 




Figure 6: Higgs Penguin contributions to the processes (a) r 3fi and {b)Bs — > r/i. The 
effective vertices gts and are generated from loops involving SUSY partners which are 
then integrated out. The mediating 5'° stands for neutral Higgs mass eigenstates, h^, H^, A^. 



Going one step further, the lepton flavour violating Higgs couplings can also be combined 
with the quark flavour changing coupling studied earlier. In this way we can also study the 
LFV and FCNC process r^. In the MSSM with large tan/3 the dominant Higgs 

contribution will again come from the penguin diagram mediated by the Higgs bosons as 
shown in fig. 6b. The branching ratio for this decay may be written as. 



Bt(B. 



2567r 



f2 

,4 JB, 



2 

mt 



1 2 



Ml, 



+ 



+ 



X 



(1 - a;^ - Xr)\Jl- 2{xr + x^) + {xr - Xf,y. 



(5^ 



Here we have defined x,- 



[mij Mb^Y ■ We have concentrated here on the final state ^ 

2 „2 l2 



because the contribution for the final state r yU+ goes like, -^^^ + -jr^ — , and approx- 

imately vanishes at large tan p. It is worth noting that at present there are no experimental 
bounds set on the process Bg —>■ rfi and only weak bounds set on the related processes. 



BR(5rf 
BR{Bd 
BR(5, 



fie) 

Tfl) 

■ re) 



< 
< 
< 



1.8 X 10"^ at 90% C.L. 
3.8 X 10"^ at 90% C.L. [34] 
3.8 X 10"^ at 90% C.L. [34] 



In the following section we shall discuss the numerical procedure used in our analysis 
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of the above phenomenology. The full list of input parameters and electroweak observables 
used in the numerical fitting are also given in detail. 

4 Numerical Analysis 

In our numerical analysis we have adopted a complete top-down approach [35]. At the GUT 
scale the MSSM gauge couplings are related to the GUT scale couplings as = = acuT 
and QJs — Q;G[/r(l + 63), where €3 sums up the effects of GUT scale threshold corrections. We 
shall be studying the SU (5) theory as outlined in section 2. The complete list of parameters 
of the Yukawa sector at the GUT scale are, 

/ diX^ \ / uiX^ \ 

= d^X' , Yf^^ = yj^^ u,X' VcKM, (59) 
\ ds J \ Us J 

( mX^ \ / A^ilO'^ \ 

1;^'"" = C/pMNS n^X^ , Mj^^= iVslQii .(60) 
\ 0.6 / \ N^W^ I 

Here, A = 0.22, and the coefficients, di,ei,Ui,ni and Ni are all of order 1 and are varied 
during our analysis in order to obtain correct fits to electro-weak data. We have chosen to 
keep Y^^^^ fixed at 0.6 in order to ensure interesting rates for r — > fi'j. In fact, in 5*0(10) 
unification a neutrino Yukawa coupling of this size would certainly be expected. In addition 
the three CKM and PMNS mixings are defined as follows, 

sin = 0.22 ai2, sin ef^ = 0.0026 «i3, sin 9^^ = 0.03 aas (61) 
sinei^^^Pu: sin^i3-0, sin^^3-^/?23. (62) 

Again the parameters ctjj and Pij are of order 1 and through their variation a good fit to the 
low-energy CKM and PMNS matrix elements can be obtained. 

A complete list of input parameters of our model are then as follows, 

MguTi OiGUT, £3, dl,2,3, '"1,2,3, ^1,2, -^1,2,3, Q;i2, 13,23, (^12,23 

(63) 

Ml/2, ^0, KMz), tan/3, mo, and m^. 

The top-down approach implies that we can freely vary or hold fixed any one of them and 
then investigate the fit properties. This is one of the advantages of doing the analysis top- 
down. 

We note that taking advantage of the top-down approach we have kept Aq — 0, ii{Mz) — 
120 GcV and tan/3 = 50 fixed throughout the analysis. Here we are taking the hint from 
the Brookhaven muon g — 2 experiment that the sign of should be positive [36]. 

Two-loop RGEs for the dimensionless couplings and one-loop RGEs for the dimensionful 
couplings were used to run all couplings down to the scale Al^f> where the heaviest right- 
handed neutrino is decoupled from the RGEs. Similar steps were taken for the lighter M2R 
and Miji scales, and finally with all three right-handed neutrinos decoupled the solutions 
for the MSSM couphngs and spectra were computed at the Z scale. This includes full one 
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Observable 


Experimental 


a 


Name 


Central value 






137.04 


0.14 




0.11664 X 10"" 


0.00012 X 10"" 




0.1180 


0.0020 


Mt 


172.7 GeV 


5.1 GeV 




1.25 GeV 


0.20 GeV 


rriii 


2.75 McV 


1.25 MoV 


nil, 


4.250 GeV 


0.20 GcV 


nis 


0.105 GeV 


0.035 GcV 


rud 


6.000 MeV 


2.0 MeV 


Mr 


1.7770 GeV 


0.0018 GeV 




0.10566 GeV 


0.00011 GeV 




0.51100 MeV 


0.00051 McV 


Vu, 


0.2257 


0.0022 


Vcb 


0.0413 


0.0015 


Vub 


0.00367 


0.00047 


Mz 


91.188 GeV 


0.091 GcV 


Mw 


80.419 GeV 


0.080 GeV 


Pnew 


-0.2 X 10-3 


0.011 X 10-3 


Br(6 ^ 57) 


0.3470 X lO--'' 


0.045 X 10-3 


Afi,, 


0.245 X K)-*^ 


0.090 X 10-** 



Amit^ 0.24 X 10-2 eV^ 0.030 x 10-^ eV'^ 

sin^e Atm 0.4400 0.018 

Srnl^i 0.792 x 10-* eV"^ 0.035 x lO"" eV^ 

sin^ 8g„i 0.314 0.029 



Tabic 1: Tabic of obscrvablcs and a values which are used to calculate the that enables 
best fit regions to be determined via minimisation. Masses denoted with a capital M are 
pole masses while a lower case is used for MS running masses. The bottom quarks running 
mass is evaluated the scale Mf, whilst the light quarks are evaluated at the scale 1 GeV. 



loop SUSY threshold corrections to the fermion mass matrices and all Higgs masses while the 
sparticle masses are obtained at tree level. A detailed description of the numerical procedure 
can also be found in reference [35]. 

In our approach and were varied to optimise radiative elcctroweak symmetry 
breaking (REWSB), which was checked at one loop with the leading and mf corrections 
included following the effective potential method in [37] . 

As mentioned earher, the calculation of the MSSM Higgs widths are performed with the 
numerical code HDECAY[30]. This FORTRAN code allows Higgs widths and branching 
ratios to be calculated with QCD corrections for hadronic channels included. Here we make 
use of the total width of the neutral Higgs to enable the branching ratios for lepton flavour 
violating neutral Higgs decays to be evaluated. 

An experimental lower bound on each sparticle mass was imposed. In particular, the most 
constraining are: the LEP limits on the charged SUSY masses (m^±,mf > 105 GeV), the 
CDF limit on the mass of the CP odd Higgs state (m^o > 105-110 GeV, vahd for tan/3 
50) [38], and the requirement that the lightest SUSY particle should be neutral. Finally, 
a function, {Xf'^ — X^^^)'^/a^ , is evaluated based on the agreement between the 
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400 600 800 1000 1200 1400 



mo [GeV] 

Figure 7: Contour plot of BR(r fi^y) in the ttiq — M1/2 plane with /i = 120 GeV and 
Ao = 0. Contours are in units of 10"'', the present experimental bound marked as "3.1". 



theoretical predictions and 24 experimental observables collected in table 1. The numerical 
procedure described above was previously exploited to study lepton flavour violation, the 
muon anomalous magnetic moment, neutrino oscillations and 6 — >^ 57 within supersymmetric 
Pati-Salam[21, 39] and 50(10) [35, 40] models. 




Figure 8: Variation of the BR(r 
ranges from 0.1 to 1.0. 



/i7) as the third generation Neutrino Yukawa couphng 



For our analysis we selected 100 evenly spaced points in the SUSY parameter space each 
with values of mo and M1/2 in the plane, mo = 300 - 1400 GeV and M1/2 = 300 - 1000 GeV. 
For each of the 100 points we also held fixed /i = 120 GeV, tan/? = 50 and Aq = 0. The 
remaining input parameters of eq. (63) are allowed to vary in order to find the minimum 
of our electroweak function. The results of our analysis are discussed in the following 
section. 
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5 Discussion 



In the previous section we outlined the numerical procedure used to analyse 100 points in 
the SUSY mo — M1/2 plane. In our analysis of these points we have kept fixed // = 120 GeV, 
tan P — 50 and Aq — 0. In this section we shall present in detail the results of our analysis. 



5.1 Lepton flavour violating /j, ^ ej and t ^ /ij 

Let us firstly look at the important lepton fiavour violating decay t fi'~f. A contour plot 
showing the variation of the branching ratio for r — > /X7 in the mo — M1/2 plane is given 
in fig. 7. This plot shows that the rate for this process is very high over most of the plane 
and even in the low mo — M1/2 region the rate can exceed the present experimental upper 
bound. These large rates are a direct consequence of our choice of Y„.^ = 0.6. It is interesting 
therefore to notice the variation of the rate for r ^ fi^y with the size of Fj,.^ as shown in 
fig. 8. In this plot we have allowed Y^^ to vary from 0.1 to 1.0 whilst observing the resulting 
change in the rate for t fj.'j. We observed this variation for a small number of fixed points 
in the mo — M1/2 plane, as noted in the figure. It is interesting to notice that the decay 
rate changes by about three orders of magnitude whilst the Yukawa coupling varies by just 
a single order of magnitude. This type of Yi, dependence clearly comes from the dependence 
of Am| shown eq. 24, which shows that (Am|)23 has a roughly (Y^^^^)"^ dependence. 
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Figure 9: BR(/x —>■ ej) plotted against the Pseudoscalar Higgs mass, m^o. 



On the other hand, the rate for the decay ji ^ appears to be much more suppressed. 
Our results, plotted in fig. 9 against the Pseudoscalar Higgs mass, show that for small values 
of m^o the rate is highly variable over the range 10^^^ - 10^^^ At larger values of 
a slight bunching of points appears to occur in a narrow horizontal band. It seems that 
our model prefers a ^ ^ 67 rate in the region lO""*^^ — 10~^^, which is close to the search 
limitations of the present MEG experiment [8]. 



5.2 Higgs mediated Bg — ^ Bg — ^ r//, t ^ S/j, and AMg 

The plots for the Higgs mediated contribution to the fiavour changing neutral current decay, 
Bs — > jifJ,, are presented in fig. 10b. As previously reported by a number of authors[13, 26] 
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the branching ratio for Bg /i/i is particularly interesting with rates ranging from 10~^° for 
heavy m^o up to almost 10~^ for a particularly light m^o. The present 90% CL. experimental 
upper bound of 0.8 x 10^^[14] excludes the highest predicted rates and hence is beginning 
to probe the Higgs sector into the region rriAo ~ 150 GeV. This rare decay is certainly of 
particular interest and future studies will continue to probe the Higgs sector to an even 
greater extent. Within the standard model the expectation is that Bt{Bs —>■ fJ^fi)sM ~ 10~^. 
Therefore we can see from fig. 10b that the Higgs mediated contribution to this process 
would dominate as long as m^o < 500 GeV and this decay has the possibility of being the 
first indirect signal of supersymmetry. 
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Figure 10: The Higgs mediated contribution to Bt{Bs — > /i/i, r/i) plotted against the Pseu- 
doscalar Higgs mass. In figure (b) the horizontal line at 0.8 x 10~^ shows the present 90% 
CL. experimental upper bound [14]. 



The lepton flavour violating decay Bg —>■ Tjj, is plotted in fig. 10a against the Pseudoscalar 
Higgs mass. This plot shows that the branching ratio for this decay could be as large as 
~ 10~^°, but is certainly not as interesting as the flavour conserving decay just discussed. 
The branching ratio for the decay r — /i/i/i is of a similar order of magnitude to that in 
fig. 10a. These rates are not quite as high as those reported in eq. (12) and (13). 

In the limit of large tan /?, Bg fi~^^~ and AMg are correlated. This correlation is shown 
in the two panels of fig. 11. For these two panels the two different values of listed in 
eq. (45) and (46) are used. The upper panel (/^^ = 230 MeV) shows that the central value 
of {AM^^^ - AMf^) coincides with the bound from Br(5^ IX^ jjL ). The lower panel 
[Jbs = 259 MeV) shows that the data points with AMf ^ at the central value, are in fact 
ruled out by the bound on Br^Bg —> n'^fj,~). The uncertainty in the SM prediction for AMg 
is rather large and in fact all of the data points of fig. 11 agree with the recent Tevatron 
measurement to within la. These two panels clearly show that the interpretation of the 
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recent measurement depends crucially on the uncertainty in the determination of Jb. 



Figure 11: The correlation of Bt{Bs ^^A*^) and AM^ are given for /b^ = 230 MeV 
(upper panel) and /b^ = 259 MeV(lower panel). The horizontal line in both panels shows 
the present experimental bound on Br(i?s — * M^^~)- The vertical lines mark the central 
value of (AMpDF - AMf^O- 



Figure 12: Plot for the ratio, AM^ / AM^, of the Higgs contribution to the operators 
Qf^^ and Qf^. 



The plot in fig. 12 shows the ratio, AM^^ /AM^^, of the contributions to the operators 
Qf^^ and as defined in eq. (44). It is commonly assumed that the contribution to the 
qSll operator, AM^^, is negligible. From fig. 12 we can see that the contribution of the 
qSll operator, AM^^, is between 40% and 90% of AM^^ and hence is significant. 



5.3 Lepton flavour violating decays of MSSM Higgs bosons 

Our numerical results for the branching ratio of r/x are presented in fig. 13c. The 

rates for this process are plotted against the Pseudoscalar Higgs mass, m^o. The predicted 
rates range from 10~^ up to a few times 10"''. In fig. 13c we can see a broad band of 
points stretching from m^o = 100 GeV up to m^o = 600 GeV showing that most points 
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are predicting a rate of around 10^^. The decay rate appears to be almost independent of 
the Pseudoscalar Higgs mass, although there is a slight peak at the lower range of the Higgs 
mass, which is where the highest rates are achieved. Notice that the data points shown in 
fig. 13 are divided into two groups. The grouping is determined by whether the Bg fj^fi 
bound is in excess or not, as indicated in the figure caption. From this grouping we can 
see that the points with the highest predicted — > r/i rates appear to be excluded by the 
Bg — > fifj, bound. In this way the Higgs mediated contribution to the Bg decay can provide 
additional information on the allowed Higgs decay rate. 
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Figure 13: BR{h'^ , A'^ , r/i) plotted against ttiao- Each of the scatter points represent 

the output result from an fixed value in M1/2 and mo. For each of these fitted points we 
have assumed, /i = 120 GeV and Aq = 0. Here the dots(diamonds) mark points for which 
Br{Bs nn) is below(above) the present experimental limit. 



Fig. 13b shows our predictions for the lepton flavour violating Pseudoscalar Higgs decay, 
A^ —>■ Tfi. The rates for the decay of the Pseudoscalar are almost identical to those of the 
heavy CP-even Higgs shown in fig. 13. 

The rate for the decay of the lightest Higgs boson is shown in fig. 13a where it is again 
plotted against the Pseudoscalar Higgs mass. The predicted rates for this decay show a 
very different dependence upon m^o, as they are spread over a large range from 10^^^ to 
10"''. The branching ratio for the lightest Higgs appears to be inversely proportional to 
the Pseudoscalar Higgs mass. Hence this LFV decay will only be interesting if m^o < 300 
GeV where its rate can be comparable to those for the other neutral Higgs states. The 
data points for the lightest Higgs boson decay of fig. 13a have again been divided into two 
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groups. As before, the two groupings depend upon whether the Bg — * /i/i bound is being 
exceeded or not. The plot shows that the data points for particularly light Pseudoscalar 
Higgs mass are excluded by the Bs — > yU/i bound. These excluded points correspond to the 
largest predictions for the decay Tfi and leaves, 3 x 10~^, the highest allow decay rate. 
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Figure 14: Br(if° rfi) plotted against Br( 



Fig. 14 shows the correlation between the two LFV decays — > r/i and r — > fi'y. The 
plot shows a very nice linear relation between the two branching ratios. The rate for r — /X7 
is generally a factor of ~ 10 larger than that for Tfi, although there is a region towards 

the bottom left of the plot where the two rates become comparable. 
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Figure 15: Bt{H Tfi) and Br(r —>■ fi'-f) plotted along the line mo = M; 



1/2- 



In some other works [12] a non-decoupling effect in the Higgs decay ^ Tfi was ob- 
served. It was found that the decay t ^ fij experiences a decoupling suppression for large 
values of mo and Mi/2- This decoupling is experienced as a result of heavy SUSY particles 
in the loops of lepton flavour violating diagrams. On the other hand the Higgs LFV coupling 
do not experience such a large suppression. It was noted that, in this decoupling region, the 
LFV Higgs decay rate can be larger than the tau decay rate. 

In order to examine such an effect we present a plot of the two decay rates against 
mo = Ml/2, see fig. 14. We see that the rate for the tau decay is initially an order of 
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magnitude greater than that for the Higgs decay. As we move from small to large values of 
mo = Ml/2 the rate for the tau decay falls sharply where as the rate for the Higgs decay falls 
less sharply. By the time we reach mo = M1/2 — 2000 GeV we find that the two rates have 

become comparable and the Higgs decay rate may even become marginally larger than that 
for the tau decay. Other authors noted a greater shift than we have found, with the Higgs 
decay rate becoming ~ 6 times greater than the tau decay rate. Here we did not see such 
an large effect at large values of mo = M1/2. 



6 Conclusions 

We have analysed a supersymmctric model constrained by SU{^) unification with right- 
handed neutrinos providing light neutrino masses via the seesaw mechanism. We have been 
concerned with making predictions for lepton flavour violating decay processes such as r — > 
/i7, // — > 67 and 0° — > r/i. With lepton flavour violation in mind, we chose to study the 
large tan (3 region of parameter space where such effects are enhanced. In addition we chose 
yj,., ~ 1 under the influence of 5*0(10) unification in which tan/? is naturally large. 

Our numerical procedure utilises a complete top-down global fit to 24 low-energy 
observables. Through this electroweak fit we were able to analyse lepton flavour violating 
decay rates for rare charged lepton decays. The choice of a large third family neutrino 
Yukawa coupling naturally leads to large 23 lepton flavour violation. Within this scenario 
we have also used our numerical fits to analyse the branching ratios for the FCNC processes 
Bg nn, Bg TyU and the lepton fiavour violating decays of MSSM Higgs bosons. As a 
result, we have been able to make realistic predictions for these decays while ensuring the 
maximum r — /i violation allowed by the present r /i7 constraint. 

Our model predicts a r — > /X7 rate in the region (10~^ — 10~^) and a /x — > 67 rate in 
the region (10~^^ — 10~^^). We have seen that the branching ratio for — > Tji could 
be particularly interesting with a rate as high as few 10~^. We also saw that in the large 
mo — Ml/2 region of parameter space this branching ratio could become comparable to that 
for r ^7. The Higgs mediated contribution to the branching ratio for Bs was 
found to be particularly large and could even exceed the current experimental bound. In 
addition to the constraint of r /X7, we also saw that the Bg fj^fj, bound may also act as 
a stringent restriction on the allowed rates for 0° — > t/j,. We found the rates for the LFV 
decays Bg t/i and r — > /i/iii are rather small < 10~^°. 

The correlation of Bg — > and AMj in the large tan P hmit was also studied. The 

constraint from the recent Tevatron measurement is highly dependent upon the determina- 
tion of Jb^- It was found that the Higgs contribution to the often neglected operator Qf^^ 
may be as large as 90% of the contribution from the operator Q2^. Therefore the operator 
Qf^^ is certainly non-negligible. 
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